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基于 双向 流 固 耦合 方法 的 东风 局 特性 数值 研究 
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摘 EO 风扇 转子 工作 在 高 压 比 、 高 转速 、 高 负荷 的 级 环境 中 ,实际 流动 中 存在 风扇 的 长 叶片 与 流体 间 强 烈 的 耦合 作用 . 
针对 风扇 实际 工作 中 存在 的 流 固 耦 合 问题 ， 本 文通 过 基于 时 域 推进 的 双向 流 固 磷 合 方法 对 风扇 转子 的 流 固 耘 合 问题 进行 
求解 , 并 在 此 基础 上 对 该 转子 的 超声 速 非 失速 颜 振 现象 进行 研究 . 分 析 结果 表明 , 考虑 灯 合 效应 的 转子 特性 线 整体 向 大 流 
量 方向 偏 移 0.7%， 效 率 降低 约 0.76%， 总 压 比 基本 保持 不 变 ， 稳定 工作 裕 度 略 有 减 小 ; 转子 在 设计 转速 下 近 效 率 最 大 点 
处 , 可 能 会 发 生 基于 叶片 一 阶 固有 模 态 振动 的 超声 速 非 失速 额 振 现象 . 
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Numerical Investigations on Characteristic of Fan Based on Two-Way 


Fluid-Structure Interaction Approach 


ZHANG Shuail2 GAO Li-Min'? ZHENG Tian-Longl? DENG Wei-Min!? 
(1. School of Power and Energy, Northwestern Polytechnical University, Xi'an 710072, China; 
2. Collaborative Innovation Center of Advanced Aero-Engine, Beijing 100191, China) 
Abstract There is strong coupling between the long fan blades and the fluid in actual flow field, 
because the fan rotor works in the high pressure ratio, high rotation speed and high load level 
environment. In this paper, the fluid-structure interaction(FSI) problem of the fan rotor is solved 
by two-way fluid-structure interaction method based on the time marching, and the supersonic 
unstalled flutter phenomenon of the fan rotor is studied. The numerical results show that the fan 
rotor characteristic line shift 0.796 towards the large flow direction, the efficiency is reduced 0.76%, 
the total pressure ratio remains the same and the stable operation margin decreases slightly when the 
coupling influence is considered. The supersonic unstalled flutter based on the first-order natural mode 
vibration of the blade may occur near the maximum efficiency point at the design speed of the fan rotor. 
Key words fan; fluid-structure interaction; geometric deformation; aerodynamic performance; un- 
stalled flutter 


0 引 言 


随 着 现代 先进 航空 动力 装置 性 能 的 迅猛 提升 ， 。” 题 的 研究 大 多 是 通过 流 固 耦 合 的 分 析 方 法 来 实现 的 ， 
航空 发 动机 向 着 高 推 重 比 、 高 负荷 、 高 性 能 的 方向 。” 即 采用 耦合 求解 非 定 常 气动 力 和 结构 体 瞬 时 变形 的 
不 断 发 展 。 对 于 承担 工 质 加 功 功能 的 风扇 /压气 机 ， ”方法 ,这 种 方法 能 够 提高 数值 模拟 的 精度 ， 特 别 是 
级 压 比 越 来 越 高 ， 叶 片 长 度 越 来 越 长 ， 厚 度 越 来 越 。 在 非 线性 影响 强度 较 强 的 跨 声速 区 域 ， 因 为 该 方法 


jl 


薄 ， 叶 尖 切 线 速 度 越 来 越 高 ， 非 定常 气动 负荷 不 断 能够 反映 流体 流动 与 结构 振动 的 耦合 物理 过 程 ， 不 
增 大 .叶片 在 气动 力 和 离心 力 共 同 作用 下 的 流动 不 ”但 考虑 到 结构 变形 引起 气动 力 变化 的 非 线性 ， 还 可 
稳定 问题 、 气 动弹 性 问题 等 更 加 突出 ;而 叶片 变形 ” 以 处 理 结构 动力 学 方面 的 结构 体 大 变形 、 几 何 及 材 
后 使 得 叶 型 偏离 最 佳 理 想 设 计 而 影响 风扇 /压气 机 。” 料 非 线 性 等 问题 O. 而 随 着 计算 机 技术 的 迅猛 提升 ， 


的 气动 性 能 ， 严 重 时 会 对 叶片 的 结构 强度 及 振动 计算 流体 力学 /计算 固体 力学 (CFD/CSD) 耦合 计算 
性 产生 较 大 影响 了 方法 也 得 到 了 迅速 的 发 展 B, Wilson 等 9 用 非 线 

目前 国 内 外 针对 叶轮 机 械 中 的 叶片 气动 弹性 问 性 气动 弹性 模型 研究 了 跨 声速 风扇 叶片 在 不 同 工 作 
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状态 下 ， 几 何 变形 对 性 能 的 影响 及 叶 尖 相对 安装 角 ”求解 方法 、 油 流 模型 等 的 正确 及 合理 性 ， 对 该 风月 
的 变化 , 发现 叶 片 在 不 同 转速 下 的 变形 对 气动 性 能 ”转子 的 整 机 在 设计 转速 下 进行 定常 流 场 特 性 数值 模 
影响 显著 . W. 图 2 所 示 为 该 风扇 两 级 叶片 的 流 道 及 计算 网 格 . 
随 着 计算 力学 的 发 展 ， 数 值 计 算 方法 在 叶轮 机 ”风扇 的 单 通道 四 排 叶 片 网 格 总 数 150 万 ,不 考虑 各 
械 叶片 设计 过 程 中 扮演 着 重要 的 角色 C 可 以 有 效 。 级 转 、 静 子 的 叶 顶 间隙 ， 其 进出 口 条 件 保 持 和 实验 
预测 叶片 在 实际 工作 过 程 中 存在 的 颜 振 等 气动 弹性 ”环境 条 件 基本 一 致 ,两 级 转子 的 转速 均 为 设计 转速 ， 
问题 S, 本 文采 用 基于 时 域 推进 的 双向 流 固 耦 合 数 ”进口 总 压 为 101325 Pa, 总 温 为 288.15 K, 同时 对 第 
值 计算 方法 对 NASA CR-120859 风扇 第 一 级 风扇 转 ”二 级 静 子 的 背 压 进行 调节 ; 气体 动力 学 方程 的 求解 
子 在 设计 转速 下 的 气动 弹性 问题 进行 数值 计算 研究 ， 采用 求解 定常 的 雷诺 时 均 Naiver-Stokes(N-S) 方程 
探究 风扇 转子 叶片 的 静 气 动弹 性 变形 规律 及 其 对 风 的 方法 ,清流 模型 采用 考虑 气体 前 切 应 力 的 SST i 
扁 转 子 性 能 、 内 流 场 的 影响 ; 并 在 此 基础 上 , 进一步 。 流 模型 。 
研究 该 风 遍 转 子 的 超声 速 非 失速 颤 振 现 象 ， 为 叶片 
设计 和 风扇 实际 工作 中 避免 此 类 气动 弹性 不 稳定 现 
象 提供 数值 计算 与 理论 借鉴 


1 计算 模型 建立 


1.1 算 例 验证 
本 文 所 采用 的 数值 模拟 理论 及 方法 均 基 于 


NASA CR-120859 风扇 ， 此 风扇 是 包括 适度 的 高 

RER, RADR, SAIRAS (00 H2 AEAN 

计 概 念 的 两 级 风 B. 该 风 转子 的 第 一 级 转子 叶片 Fig. 2 The calculation grid of two stage fan 

的 弦 长 较 宽 、 展 弦 比 较 大 、 叶 片 较 薄 并 且 带 防 振 吓 

站 第 二 级 转子 的 叶片 较 类 入 度 较 大 、 展 纱 比 相对 PHA S 
较 小 , 该 风扇 的 各 级 动 、 更 叶片 的 实物 如 图 1 所 示 ， 对 比如 图 3、 图 4 所 示 . 从 对 比 图 中 发 现 , 实验 测量 


与 数值 计算 的 绝热 效率 与 总 压 比 在 变化 趋势 上 完全 
吻合 ， 在 数值 上 略 有 差异 ， 实 验 测量 的 绝热 效率 略 
高 于 数值 计算 值 ， 而 数值 计算 的 总 压 比值 略 高 于 实 
验 测量 值 . 实验 测量 的 出 口 测 点 位 于 第 二 级 静 子 出 
口 截面 ， 而 数值 计算 的 效率 计算 截面 位 于 风扇 整 机 
的 延伸 段 出 口 ， 这 导致 实验 测 基 的 静 温 值 大 于 数值 
计算 ，, 实验 测量 绝热 效率 略 大 于 数值 计算 值 , 在 几 
可 建 模 过 程 中 对 该 风扇 进行 简化 ， 数 值 计算 过 程 中 
不 考虑 风扇 转子 的 第 一 级 转子 的 叶片 防 振 突 肩 (如 


a 
T 
C 


(a) First-stage blade (b) First-stage vane 图 1(a) 所 示 ) 及 各 级 转 、 静 子 的 叶 顶 间隙 等 因素 , 这 
些 均 可 能 使 该 风扇 转子 的 数值 计算 总 压 比 大 于 实验 
WEA. 


综 上 所 述 ， 本 文 所 选取 的 数值 计算 方法 、 满 流 
模型 、 转 子 几何 模型 简化 等 经 过 检验 , 风扇 转子 的 气 
动 性 能 和 实验 测量 的 变化 趋势 及 数值 基本 吻合 ， 能 
够 满足 本 文 对 于 转子 叶片 看 合 特性 研究 的 要 求 ， 
于 受 困 于 数值 计算 的 计算 量 、 转 静 交 界面 数据 传递 、 
动 欧 干涉 等 因素 ， 本 文 在 接 下 来 的 叶片 契合 特性 研 


(c) Second-stage blade (d) Second-stage vane 

M T 人 Hbk HA o B 
图 1 NASA CR-120859 RRUA ETAIT LYE 究 中 均 选 用 考虑 叶 顶 间隙 、 不 带 减 振 突 肩 的 NASA 
Fig. 1 Blades physical pictures of NASA CR-120859 CR-120859 风扇 的 第 一 级 转子 . 


1.2 本 文 计算 模型 
为 验证 本 文 所 选取 的 几何 模型 简化 方法 、 流 场 NASA CR-120859 风扇 第 一 级 转子 叶片 长 达 
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8 期 


0.236 


dl. 


TJ 


实际 转动 过 程 中 


响 变形 大 、 


m， 进 口 轮 载 比 达 到 


比 风 


P, KA 
耦合 特性 明显 ， 


张 帅 等 : 基于 双向 流 固 


.4， 叶 片 的 设计 参数 见 


户 第 一 级 转子 叶片 展 弦 比 大 , 在 转子 


Ei 叶 


的 基础 上 而 


特性 的 影响 


究 静 气动 变 
， 选 择 该 两 级 风 


究 对 象 , 该 转子 的 三 维 几何 模型 如 图 5. 


片 受 离心 力 及 气动 力 影 
因此 为 在 流 固 耦合 分 析 
形 对 风扇 性 能 、 流 场 及 振动 
扇 的 第 一 级 转子 作为 研 
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Fig. 3 Comparison of the fan adiabatic efficiency 
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4 风扇 整 机 总 压 比 对 比 图 


EET f 


1.01 


Fig. 4 Comparison of the fan total pressure ratio 


zi 


格 生成 是 数值 计算 中 很 重要 的 环节 ， 其 质量 


决定 着 数值 计算 能 否 完成 及 其 结果 是 否 可 靠 ， 而 双 


TET RE D 


耦合 数值 计算 对 于 流体 域 及 固体 域 的 网 格 质 


合 方法 的 某 风扇 特性 数值 研究 


推进 的 非 定 常 计算 ， 计 算 域 的 网 格 是 
变化 的 动 网 格 ， 对 于 本 文 而 
转子 叶片 ， 其 叶 尖 变 
杂 ， 在 数值 计算 过 程 中 很 容易 出 现 负 
值 计算 无 法 继续 进行 。 经 过 大 量 的 网 
数量 级 网 格 的 数值 计算 发 现 , 流体 域 
下 比 较 敏 感 ， 在 数值 计算 过 程 
而 非 结构 化 网 格 的 抗 畸变 能 力 相 


变 
格 ， 


究 的 带 叶 
E RB, 


随 着 时 间 不 断 
硕 间 院 的 风 
叶 顶 间隙 网 格 复 


TR 3 


格 ， 导 致 数 
型 及 不 同 
网 格 对 于 叶片 
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Pa 


PIA 
对 


易 产 生 负 网 
EER, XIF 


风扇 在 近 失 速 或 是 近 堵 塞 的 工 况 都 能 很 好 的 适应 ， 
流体 域 的 非 结 构 化 网 格 节点 数 为 621354， 单 元 数 为 


2322831， 固 壁 处 附 下 


层 层 数 为 15, W 


体 壁 


层 网 格 尺度 满足 求解 方程 对 于 v^ 值 的 要 求 。 
表 1 NASA CR-120859 第 一 级 风 扁 转子 主要 参数 


Table 1 The main parameters of the first stage 


面 的 第 一 


fan rotor 

设计 参数 数值 

叶片 数 28 
叶 高 /m 0.2361 
设计 转速 /(r/min) 10720 
叶 尖 切线 速度 /(m/s) 441.96 

叶 顶 间隙 /mm 1.5 
设计 流量 /(kg/s) 83.55 
设计 效率 0.894 
设计 压 比 1.787 


图 


本 文 基 于 


寺 域 推进 的 双向 流 固 
解 非 定常 雷诺 平均 N-S 方程 
体 ， 瞬 态 求解 格式 是 基于 动 
后 差分 格式 , 详细 求解 器 参数 设置 如 表 2, 


5 第 一 级 转子 三 维 几 何 模型 


Fig. 5 The three dimensional geometric model of the first 


stage rotor 


IR 


上 


固体 


的 风 


基 要 求 更 高 .双向 流 固 耦 合 的 计算 方法 是 基于 时 域 


解 基于 瞬 态 动力 学 方程 


域 的 结构 瞬 态 分 析 采 


形 对 流体 域 流 场 的 影响 ， 


性 主要 关注 的 是 四 


耦合 解法 采 
， 流 体 介 质 选 取 理 想 气 
网 格 方法 的 二 阶 欧 拉 向 


男 体 


体 


因此 选取 先 求解 流体 


IR 


BK 


的 直接 积分 方法 . 在 双 
的 耦合 时 域 推进 方法 中 ， 流 体 域 的 非 定常 计算 与 


j 相 同 的 时 间 步 长 ， 而 流 
体 域 与 固体 域 的 瞬 态 求解 是 顺序 进行 的 . 本 文 研 
扇 转 子叶 片 的 耦合 特 ; 


究 
BE 


成 再 
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求解 固体 域 的 顺序 
表 2 求解 器 参数 设置 


Table 2 The parameter settings of solver 


计算 参数 数值 
济 流 模型 SST 
最 小 时 间 步 (t/s) 1,0x10— 
总 时 间 (T/s) 0.12 
虚拟 时 间 步 (Step) 10 
参考 压力 (Pief/Pa) 101325 
进口 总 温 (Tr1/K) 288.15 
进口 总 压 (Pri/Pa) 0 
出 口 背 压 (Ps /Pa) Variables 
2 风扇 转子 耦合 特性 分 析 
对 NASA CR-120859 风扇 转子 进行 非 定常 计算 
和 基于 时 域 推进 的 双向 流 固 耦合 计算 ， 对 计算 得 到 


的 转子 性 能 参数 进行 时 均 化 处 理 , 图 6、 图 7 所 示 为 
两 者 的 总 性 能 参数 对 比 图 . 叶片 的 气动 弹性 变形 导 
致 转子 的 等 转速 线 发 生 明 显 的 偏 移 ， 耦 合 后 的 特性 
1 线 整体 向 大 流量 方向 偏 移 0.7% 左 右 ， 变 化 趋势 旨 
本 吻合 . 效率 数值 有 所 降低 , 最 大 效率 值 降低 0.76%, 
耦合 前 后 的 总 压 比 在 数值 上 基本 相同 ， 转 子 的 稳定 
I 


E 裕 度 略 有 减 小 . 在 双向 流 固 耦合 计算 过 程 中 ， 
于 随 着 时 间 的 推进 , 叶片 发 生 弯 扭 振动 , 导致 转子 叶 


片 的 各 截面 安装 角 变 小 , 进口 气流 角 变 大 , 流 道 的 通 
流 面积 发 生变 化 ， 使 得 风扇 转子 的 特性 曲线 发 生 整 
体 偏 移 ; 同时 , 由 于 转子 叶片 的 往复 振动 、 流 场 的 不 
断 变化 及 叶 项 间 隙 的 变化 等 因素 使 得 风扇 转子 的 等 
炉 效 率 相 比 于 耦合 之 前 有 所 降低 . 

0.94 
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Fig. 6 The entropy efficiency of fan rotor between the 


unsteady and coupling way 
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到 7 耦合 前 后 风扇 转 子 总 压 比 


Fig. 7 The total pressure ratio of fan rotor between the 


unsteady and coupling way 


2.1 流 场 分 析 

为 研究 风扇 转子 耦合 前 后 效率 及 压 比 的 变化 ， 
对 叶片 表面 的 带 压 等 参数 进行 分 析 . 图 8 给 出 
风扇 转子 叶片 在 设计 工 况 下 ， 耦 合计 算 在 叶片 
叶 高 截面 位 置 的 表面 更 压 变化 情况 ， 其 中 黑色 实 
表示 双向 流 固 耦合 的 非 定常 计算 周期 内 静 压 波动 的 
时 均值 ， 而 虚线 和 点 划 线 表示 对 应 非 定常 
片 表面 带 压 波动 的 最 大 值 与 最 小 值 . 首先 对 比分 析 
转子 叶 表 静 压 的 时 均值 ， 在 沿 着 叶 高 方向 ， 
吸 、 压 力 面 的 静 压 包 络 线 面 积 逐 渐 增 大 ， 即 从 叶 根 
到 叶 尖 ， 转 子叶 片 的 负荷 逐渐 加 大 ， 
增强 ;， 叶 表 静 压 波 动 主要 集中 于 叶片 的 吸力 面 及 靠 


近 叶 顶部 分 ， 静 压 波动 沿 着 展 向 逐渐 增强 的 同时 ， 
从 叶片 的 前 缘 波及 尾 缘 ， 而 叶片 压力 面 仅 在 叶片 前 
缘 及 98% 叶 高 的 尾 缘 处 有 较 大 的 波动 . PET ERREUR 
压 波动 变化 及 分 布 趋势 ， 与 耦合 过 程 中 叶片 的 几何 
变形 规律 紧密 相关 ， 而 叶片 变形 主要 集中 于 叶片 的 
前 缘 及 叶片 上 半 部 分 ， 


针对 转子 较 高 截面 叶 表 更 压 波动 剧烈 的 情况 ， 
如 图 9 为 时 均 化 的 叶 表 静 压 变化 趋势 对 比 图 。 图 
9(a) 中 转子 叶片 考虑 耦合 效应 的 叶 背 静 压 突变 位 置 
( 即 激 波 位 置 ) 明显 沿 弦 向 提前 7 多 左右 ， 而 压力 突 
变 的 趋势 明显 增强 ， 考 虑 耦合 效应 的 叶 背 激 波 有 所 
增强 ， 气 流 经 过 激 波 后 的 静 压 值 基本 保持 不 变 ; 图 
9(b) 中 98% 叶 高 处 的 激 波 位 置 在 叶 背 沿 弦 向 提前 了 
10% 左 右 ， 激 波 强 度 也 有 所 增加 ， 但 激 波 后 的 更 压 
值 基本 保持 不 变 ;， 而 耦合 后 的 叶片 吸 、 压 力 面 的 静 
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面积 有 所 减 小 , 转子 叶片 做 功 全 


效率 


有 


风 


AH 


以 上 分 析 ， 考虑 双向 流 固 耦合 效 
的 激 波 强度 有 所 增加 ， 激 波 位 置 沿 纺 
前 移 ， 导 致 叶片 通道 内 流动 损失 有 所 增加 ， 局 
EH IUIS, 


所 降低 ， 


3k 帅 等 : 


基于 双向 流 固 


EHA I, 


IM TESESR 


不 


时 
转 


使 得 考虑 看 合 效应 的 风月 


固体 域 计算 结果 及 分 析 
户 转 子叶 片 在 工作 过 程 中 承受 


和 气动 力 载荷, 为 了 使 得 转子 叶片 在 耦合 过 程 中 ， 


未 收敛 到 


能 小 ， 


以 及 耦合 
历程 尽 可 能 短 ， 通 过 改变 臣 


区 敛 过 程 中 叶片 几何 变 
计算 过 程 达到 稳定 收敛 位 置 的 时 让 


子 


较 大 的 离心 力 


E Et o REUS 


体 叶片 的 初始 阻尼 比 


实现 上 述 要 求 . 图 10 所 示 为 风扇 转子 叶片 在 不 同 
中 分 析 发 现 , 不 


始 阻 


尼 比 下 结构 的 响应 历程 . 从 图 


ZH 


来 
初 
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(d) 9896 Span 
图 8 风扇 转子 在 典型 叶 高 处 表面 静 压 变化 包 线 图 


blade surface 


的 初始 阻尼 比 下 ,叶片 弹性 变形 收敛 型 


初始 阻尼 比 越 小 ， 叶 片 弹 性 变 
比 越 小 ， 
于 叶片 流体 域 网 格 质 量 要 求 越 高 ， 
越 容易 出 现 负 网 格 . 当 叶 片 的 初始 
叶片 的 静 弹 性 变形 响应 的 收敛 


形 收敛 到 稳定 幅 值 的 
时 间 历 程 越 短 ; 耦合 计算 的 收 剑 过程 中 ， 初 始 阻尼 
收敛 过 程 中 的 不 稳定 振动 


u 


敛 过 程 中 的 网 格 最 大 变形 量 远 高 于 


The static pressure change package of the typical 


相同 的 幅 值 ; 


幅 值 就 越 大 ， 大 
在 计算 过 程 中 就 
阻尼 比 为 0.1 时 ， 
最 短 , 但 其 在 收 
其 他 ， 这 对 叶片 


流体 域 的 网 格 质 量 提出 了 很 高 的 要 求 。 因此， 需要 


选择 一 个 兼顾 耦合 过 程 
变形 量 的 叶片 初始 阻 
图 11 为 初始 阻尼 


H e=0.3 


4 时 间 
2 IA Mi (e= 


叶 尖 截面 前 缘 、 中 部 及 尾 缘 (LE, Mid 


历程 与 网 格 最 大 
3). 
T, X 


记 转 子叶 
及 TE) 
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在 离心 力 载 荷 与 气动 力 载荷 作用 下 的 响应 。 从 变 ”叶片 前 缘 变 形 幅 值 来 看 ， 叶 片 在 近 失 速 的 工 况 点 变 
形 幅 值 来 分 析 ， 叶 人 尖 前 缘 的 变形 幅 值 最 大 ， 叶 尖 ERRAK. 在 转子 的 等 转速 工作 线 上 , 叶片 叶 尖 截面 
中 部 变形 幅 值 次 之 ， 叶 尖 尾 缘 的 变形 量 最 小 ， 可 最 大 变形 量 随 着 流量 的 增加 而 逐渐 减 小 ,特别 是 叶 
以 发 现 叶 片 的 振动 形式 是 扭转 ， 从 振动 的 稳定 性 来 。” 人 尖 前 缘 变 形 量变 化 更 加 剧烈 ,这 说 明 此 处 离心 力 载 
分 析 ， 叶 尖 截 面 尾 缘 最 先进 入 稳定 状况 ， 而 变形 ” 和 荷 与 气动 力 载 荷 在 叶 尖 作用 方向 相同 ， 气 动力 引起 
量 最 大 的 叶 尖 截面 前 缘 振 动 幅 值 较 大 并 且 收 敛 速 度 。 的 叶片 变形 量 随 流量 的 增加 而 逐渐 减 小 .而 在 风月 
gH. 转子 的 近 堵 塞 工 况 时 ， 离 心力 引起 的 叶片 变形 量 占 
总 变 的 比重 急剧 上 升 ， 甚 至 在 某 些 工 况 出 现 叶 

T€ - 片 的 ， F 仅 离心 力 载 荷 作 用 的 总 变形 量 ， 
uL Unsteady | 这 说 明 此 时 有 可 能 出 现 离心 力 与 气动 力 在 叶 尖 截面 
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Fig. 9 The relative pressure at the typical blade surface 
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8 其 张 帅 等 : 基于 双向 流 因 
作用 相反 的 状况 ， 气 动力 的 作用 减弱 了 离心 力 对 于 nal(Imter-Blade Phase Angle, IBPA) 变化 图 .叶片 
叶片 的 作用 ， 在 一 阶 模 态 振动 时 ， 气 动 阻尼 系数 在 IBPA=334.3° 
附近 出 现 负 值 ， 即 不 稳定 区 域 ， 而 第 二 、 三 阶 模 态 
es ————————————uua 对 应 的 气动 阻尼 系数 均 为 正 值 . 通过 计算 分 析 发 现 ， 
E A d 对 于 该 风扇 转子 叶片 , 一 阶 模 态 为 其 最 危险 的 模 态 ， 
ind: zzzi 且 其 最 危险 的 IBPA 大 概 在 334.3° 附近 . 通过 控制 
Saved IBPA=334.3° 不 变 , 改变 风 肩 转子 叶片 的 振动 频率 ， 
Bk 得 到 其 不 同 叶片 振动 频率 对 应 的 气动 阻尼 系数 ， 通 
$i [f 过 插值 得 到 转子 气动 阻尼 系数 为 零 的 非 失速 颤 振 临 
Soon. 界 点 。 
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图 12 叶 尖 监测 点 总 变形 量 随 风扇 转子 工 况 变化 


Fig. 12 The total deformation of the blade tip monitoring 
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由 文献 [9] 知 叶轮 机 械 的 主要 颤 振 形式 有 : K ri Prbew, 
速 颤 振 、 堵 寒 颤 振 、 超 声速 非 失速 颤 振 、A100 型 DO dE Uu Ra 
DR SHARA MRE, MELE FRERE - 
Rei EN eA (e K|13 /—189.34 Hz 一 阶 模 态 下 叶片 气动 阻尼 
US 中 较 具 影 响 力 的 超声 速 非 失速 颤 振 现 象 进行 Fig. 13 Aerodynamic i of blade in Pow mode 
研究 when f=189.34 Hz 


a pat dE Ac Bl Xe Te D Lp deett HH HL — 
KRAZ, KEH AEX NASA CR-120859 So 
PRAA Pe-T- Wr FrXETTAE IR MRENA. 由 [ 一 Mode2_ Acrodynamic Damping 
文献 [9], [10] 知 ,压气 机 的 非 失 速 额 振 一 般 发 生 在 转 [ aer 
子 效率 较 高 的 工 况 ， 选 取 该 风扇 转子 在 设计 转速 
的 近 效 率 最 高 点 处 的 工 况 进行 研究 . 由 叶轮 机 械 转 
子叶 片 以 往 的 实验 研究 发 现 ， 转 子叶 片 颤 振 受 其 前 
三 阶 振动 模 态 影响 较 大 ， 而 高 阶 振动 模 态 由 于 本 身 
对 应 的 叶片 振动 频率 很 高 ， 而 实际 中 叶片 发 
的 频率 较 小 且 振 动 幅 值 较 小 上 5, 所 以 文章 主要 基于 
能 量 法 52 对 该 风扇 转子 的 前 三 阶 模 态 的 振动 现象 H 
进行 着 重 分 析 . 0.003 E, 
通过 对 风扇 第 一 级 转子 的 前 三 阶 模 态 的 非 定 常 
气动 力 做 功 及 气动 阻尼 系数 计算 来 达到 分 析 转 子叶 UA Í 
片 稳定 性 问题 . 模 态 分 析 得 到 该 转子 叶片 固有 的 一 REM EF re ML 
阶 模 态 频率 为 189.34 Hz, 二 阶 模 态 频率 为 599.32 Hz, E 
三 阶 模 态 频率 为 967.69 Hz, 图 13~15 分 别 为 风 3 Fig. 14 pus itid m ed RU mode 
转子 前 三 阶 模 态 对 应 的 气动 阻尼 系数 随 叶 间 相 位 when f—599.32 Hz 
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的 形式 发 生 超声 速 非 失 速 颤 振 。 而 通过 改变 固体 
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到 15 f=967.69 Hz 三 阶 模 态 下 叶片 气动 阻尼 
Fig. 15 Aerodynamic damping of blade in third-order mode 
when f—967.69 Hz 
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Fig. 16 Aerodynamic damping of blade in first-order mode 
when IBPA=334.3° 


依据 线性 插值 对 气动 阻尼 系数 零点 进行 插值 ， 
当 气 动 阻尼 系数 为 零 时 ， 其 对 应 的 叶片 振动 频率 
为 f=201.34 Hz( 如 图 16 所 示 )， 即 为 一 阶 模 态 对 
应 的 颤 振 边 界 。 对 于 该 风扇 第 一 级 转子 而 言 ， 第 
一 阶 模 态 是 其 前 三 阶 模 态 中 最 不 稳定 的 模 态 ， 所 
以 可 以 认为 一 阶 模 态 额 振 边界 即 为 该 风扇 转子 的 非 
失速 颤 振 边界 。 以 上 分 析 表 明 ， 和 转子 在 设计 转速 
下 近 效 率 最 大 点 处 ， 一 阶 模 态 频 率 低 于 f=201.34 
Hz 时 ,转子 叶片 可 能 会 以 一 阶 模 态 、IBPA=334.3? 


叶片 


的 阻尼 (增设 减 振 突 肩 等 ) 可 以 使 叶片 避免 


非 失速 颤 振 的 发 生 ， 以 达到 使 风扇 稳定 工作 的 
Hif. 


428 de 


通过 基于 时 域 推进 的 双向 流 固 耦 合 方法 对 


NASA CR-120859 第 一 级 风扇 转子 进行 非 定常 数 
值 模拟 ， 并 在 此 基础 上 对 该 转子 的 超声 速 非 失 速 
颤 振 特 性 进行 了 研究 ,通过 计算 分 析 得 到 以 下 主要 


N 


阶 模 态 频率 低 于 f=201.34 Hz 时 ， 转 子叶 片 可 能 
会 以 一 阶 模 态 、IBPA=334.3" 的 形式 发 生 超声 速 
JEREMI. 通过 改变 叶片 的 结构 阻尼 (增设 减 
振 突 肩 等 ) 或 材料 阻尼 可 以 避免 其 非 失 速 颜 振 的 
AE. 


向 大 流量 方向 整体 偏 移 OT CAU E He uy dS 
效率 有 所 下 降 ， 最 大 效率 点 降低 0.76%; 转子 的 
总 压 比 基本 保持 不 变 ; 转子 的 稳定 工作 裕 度 有 所 


片 吸 、 压 力 面 的 静 压 包 络 线 围 成 面积 有 所 减 小 ， 叶 
片 做 功能 力 减弱 ， 


1) 考虑 双向 流 固 耦 合 效 应 的 风扇 转子 性 能 曲 


2) 考虑 耦合 效应 的 计算 结果 表明 ， 叶 背 处 激 波 


位 置 向 叶片 前 缘 移动 ， 最 大 位 移 量 达到 1096, mj 
波 的 强度 也 略 有 增强 ， 叶 片 通道 损失 增加 ; 转子 叶 


3) 转子 在 设计 转速 下 的 近 效 率 最 大 点 处 ， 一 
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